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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2669 

APPROXIMATE  THEORY  FOR  CALCULATION  OF  LIFT 
OF  BODIES ,  AFTERBODIES ,  AND  COMBINATIONS  OF  BODIES 
By  Barry  Moskowitz 


SUMMARY 

An  expression  is  developed  for  the  lift  of  a  slender  afterbody  in 
terms  of  the  "slender-body"  approximate  potential  at  the  after  end  in 
conjunction  with  a  suitably  calculated  value  of  the  potential  at  the 
forward  end.  The  failure  of  the  usual  "slender-body"  theory  to  predict 
any  lift  on  a  slender  cylindrical  afterbody  is  thereby  corrected.  The 
same  expression  is  used  to  compute  the  part  of  the  interference  lift 
generated  by  the  presence  of  a  neighboring  body,  due  to  the  interference 
upwash.  Another  expression  is  developed  to  compute  the  remainder  of 
the  interference  lift,  due  to  an  interference  pressure  gradient.  The 
lift  is  determined  for  a  cone-cylinder  body,  a  cylindrical  afterbody 
of  a  slender  wing-body  combination,  and  three  combinations' of  bodies 
to  illustrate  the  method. 


INTRODUCTION 

In  the  search  for  better  aerodynamic  configurations  for  supersonic 
flight,  the  missile  shape  has  become  more  complex.  One  of  the  problems 
arising  is  the  effect  of  interference  between  bodies  which  occurs,  for 
example,  when  the  propulsive  unit  is  mounted  external  to  the  fuselage. 

An  adaptation  of  slender-body  theory  for  obtaining  a  simple  approxima¬ 
tion  of  the  interference  lift  between  such  bodies  was  developed  at 
the  NACA  Lewis  laboratory  and  is  presented  in  this  report.  This  approxi¬ 
mate  theory  is  also  applied  for  the  calculation  of  the  lift  of  an  afterbody, 
which  is  herein  defined  as  the  portion  of  the  body  behind  the  wing  trailing 
edge  for  a  wing-body  combination  or  as  the  portion  of  the  body  aft  of  the 
nose  section  for  a  wingless  body. 

Slender-body  theory  as  originally  developed  by  Munk  in  studying 
the  lift  of  airships  (reference  l)  has  proved  useful  in  predicting  the 
lift  of  low- aspect -ratio  wings,  slender  bodies,  and  slender  wing-body 
combinations  at  supersonic  speeds  (references  2,  3,  and  4).  However, 
slender-body  theory  yields  the  unrealistic  result  that  a  cylindrical 
afterbody  of  a  wing-body  combination  and  a  cylindrical  afterbody  of 
a  body  carry  no  lift. 

Although  slender-body  theory  may  be  grossly  in  error  on  the  lift 
of  a  slender  afterbody,  the  theory  closely  approximates  the  correct 


2 


NACA  TN  2669 


» 

value  of  the  part  of  the  surface  potential  at  the  rear  of  the  afterbody  * 
proportional  to  the  angle  of  attack.  An  expression  utilizing  this  fact 
is  developed  herein  for  the  afterbody  lift  in  terms  of  the  slender-body 
potential  at  the  after  end  in  conjunction  with  a  suitably  calculated 
value  of  the  potential  at  the  forward  end.  The  more  accurate  methods  of 
calculation  of  greater  labor  (for  example,  linearized  theory)  are  thereby 
limited  to  the  forward  or  nonslender  portion  of  the  body. 

The  same  expression  is  used  to  compute  the  part  of  the  interference 
lift  of  bodies  of  revolution  due  to  the  upwash  field  generated  by  each 
body  on  the  other.  Another  interference  effect  considered  herein  is 
that  due  to  the  variable  pressure  field  generated  by  one  body  in  the 
vicinity  of  the  other,  which  results  in  a  buoyant  force. 

As  examples  of  the  method,  the  lift  of  a  cone-cylinder  body,  the 
lift  of  a  cylindrical  afterbody  of  a  slender  wing-body  combination, 
and  the  lift  of  three  combinations  of  bodies  are  determined. 


ANALYSIS 

Basic  Equations  « 

Consider  a  slender  body  of  revolution  at  zero  angle  of  attack  in 
a  uniform  supersonic  stream  of  velocity  U.  There  will  be  a  certain  * 

disturbance  velocity  field  and  an  associated  velocity  potential  produced 
by  the  body.  Because  of  symmetry,  this  velocity  field  will  give  rise 
to  no  lift. 

Now  let  the  body  be  given  a  small  angle  of  attack  a  and  be  subject, 
in  addition,  to  a  small  disturbance  field  (u,v,w)  generated  by  a  second 
body  (fig.  l).  (The  symbols  used  herein  are  defined  in  appendix  A.) 

The  first  body  will  develop  an  additional  velocity  potential  due  to  a 
and  to  the  disturbance  field.  If  v  and  w  are  suitably  restricted 
it  is  shown  in  appendix  B  that  to  a  certain  accuracy  (a)  only  that  part 
of  this  additional  potential  due  to  a  and  to  w  will  give  rise  to 
lift,  and  (b)  cross  product  terms  involving  both  potentials  may  be 
neglected  for  the  lift.  Call  this  lift-producing  incremental  potential  <P. 

Equations  of  motion  and  boundary  condition.  -  The  incremental 
potential  <j>  is  presumed  to  satisfy  the  governing  equation  for  small- 
disturbance,  isentropic  irrotational  flow 

32<PXX  -  *yy  -  ^Z  "  0  W 

In  addition,  (p  is  specified  to  satisfy  the  boundary  condition  that  the 
cross  flow,  Ua  +  w,  must  be  canceled  at  the  surface  of  the  bodyj  that  is. 
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(q>r)1^;R  *  “  (Ua  +  v)sin  9 

on  the  "body.  This  boundary  condition  is  linearized,  not  exact. 

Pressure  coefficient.  -  That  part  6p  of  the  entire  pressure 
increment  Ap  which  contributes  to  the  lift  is,  to  the  present  accuracy, 
given  by  the  linearized  Bernoulli  equation 

6p  _  -2(<PX  +  u) 

q.  u 


or 

included;  this  result  is  proved  in  appendix  B. 

Lift.  -  The  axis  of  the  slender  body  under  consideration  will  _ 
always  be  taken  parallel  to  the  x-axis.  For  simplicity  the  sidewash  v 
and  w  generated  by  the  interfering  body  are  considered  to  vary  only 
with  the  coordinate  x  in  the  general  vicinity  of  the  primary  body; 
this  restriction  is  the  one  mentioned  earlier.  (The  sidewash  will  not 
be  considered  further  because  it  contributes  nothing  to  the  lift 
(see  appendix  B)).  The  approximation  should  be  valid  if  the  bodies 
are  not  too  close  to  each  other. 


Specifically,  the  accuracy  is  such  that  all  terms  of  order  a  dR/dx 
(dR/dx)^,  but  no  higher,  which  can  contribute  to  the  lift  will  be 


(3) 


The  lift  of  the  body  may  be  obtained  by  integrating  the  component 
of  the  pressure  in  the  lift  direction  around  and  along  the  body.  The 
expression  for  the  lift  is 


where  R  is  the  radius  of  the  body  and  x^  and  X£  are  the  x-coordinates 
at  the  start  and  base  of  the  body,  respectively. 

Substituting  equation  (3)  into  equation  (4)  gives 


The  second  integral  in  equation  (5)  represents  the  lift  of  a  body 
immersed  in  a  variable  pressure  field  and  may  be  considered  as  a  buoyant 
force.  The  integral  will  be  discussed  in  detail  in  the  section  on 
interference  between  bodies  in  combination. 
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Development  based  on  first  integral  In  lift  equation.  -  This 
section  is  concerned  with  the  transformation  of  the  first  integral  in 
equation  (5)  into  a  more  useful  form  involving  the  surface  potential. 
Hereinafter  this  first  integral  will  he  called  L. 

For  constant  9,  the  quantity  Ru  dx  may  he  written 


Ru  dx  a  E  dp  -  R(<prdr)  =  d(R<p)  -  <p  ££  dx  -  R(<Pr)  ^  dx  (6) 

x  r*R  dx  rat-t'  ax 

Insertion  of  equation  (6)  into  the  first  integral  of  equation  (5)  yields 

L  *  2a  J*  sin  6 /r^  -  R^  -  ^  (7) 

0  v  x1  y 

The  equation  for  the  upwash  perturbation  velocity  is 


<pr  sin  0  +  —  cos  0 


(8) 


Substitution  of  the  value  of  <f>r  sin  6  from  this  equation  into  equa¬ 
tion  (7)  results,  upon  simplification,  in  the  following  equation  for  L: 


2it  x2  2  it 

L  -  ^  f  (*2*2  ~  Rl*l)  Sin  0  d 0  "  f  /  R§dX  f  Vd6  + 

X1 


Xo  2rt 

2a  P  8R  dx  P  S2LJ1  d.0 
U  J  dx  U  Se 


which  reduces  to 
2  it 


2  it  x2  2  rt 

L  »  y  f  (V2  -  Hlfl)  sln  0  d6  -  ^  J  R  S  ^  vf  W  d9 
0  x±  0 

because  the  integral  in  the  last  term  is  zero.  The  first  integral  in 
equation  (9)  may  be  integrated  directly  because,  as  shown  in  appendix  B, 
the  potential  is  of  the  form 

<p  a  sin  9  G(x,r) 


r 
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Equation  (9)  then  "becomes 


2q«  R2<P2  ~  Rl»l  _  p2RdEto  P 
U  sin  9  U  J  dx  J 


w  d0 


(10) 


This  equation  is  fundamental  for  subsequent  developments. 

If  and  ^2  >  the  values  of  surface  potential  at  the  forward  and 

after  ends,  respectively,  of  the  section  of  the  body  under  consideration, 
are  known  as  well  as  the  vertical  component  w  of  the  body  perturbation 
velocity,  the  lift  may  be  obtained  from  equation  (10) .  In  the  following 
sections,  procedures  are  developed  for  determining  the  lift  for  cases 
in  which  the  integration  of  the  second  term  in  equation  (10)  is  negligi¬ 
ble  or  zero.  (When  the  second  term  is  not  zero,  it  is  probably  simpler 
to  obtain  the  lift  from  equation  (5)  because  the  same  type  of  computation 
is  required  to  determine  the  w- distribution  as  is  needed  to  determine  the 
u -distribution. )  For  these  cases,  (|>2  is  assumed  to  be  the  potential  given 
by  slender-body  theory. 

Slender-body  theory  for  variable  cross  flow.  -  The  central  approxi- 
mation  of  slender-body  theory  is  the  following:  In  any  plane  x  *  con¬ 
stant,  the  value  of  the  potential  increment  <P  associated  with  the 
cross  flow  is  the  two-dimensional  potential  of  a  cylinder  having  the 
local  radius  and  moving  with  the  local  cross-stream  velocity.  For 
supersonic  flow  with  uniform  cross  flow  (that  is,  a  fixed  angle  of 
attack  and  no  variable  interference  upwash) ,  this  result  is  derived  in 
reference  3  by  taking  the  limit  of  the  linearized  value  for  radius 
approaching  zero.  The  corresponding  result  will  now  be  derived  in 
similar  fashion  for  the  more  general  case  of  variable  cross  flow. 


The  potential  in  integral  form  may  be  written 


x-Pr 

sin  0  P  g(Q(x  -  $)dg 

r  0  V(x  -  £)2  “  P2**2 


(11) 


The  expression  may  he  interpreted  as  the  potential  of  a  distribution 


of  doublets  of  strength 
appendix  B). 


f8( 


x)dx  per  unit  length 


(g(x) 


g]_ '  M  Of 


When  the  distribution  function  for  the  limiting  case  R  *♦  0  is 
related  to  the  local  cross  flow  as  in  reference  3,  there  results 


(12) 
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Inserting  this  value  of  g(x)  into  equation  (ll)  for  the  potential, 
applying  the  slender-body  approximation  (radius  approaches  zero),  and 
integrating  yield 


,  ~t  <\  R2  sin  0 

<P  -  LUa  +  W(X)J  — - - 


(13) 


On  the  hody  (r  =*  r),  equation  (13)  reduces  to 


<|>  =  Jua  +  w(x)J  R  sin  9 


(14) 


Equations  (l3)  and  (l4)  represent  the  slender-body  approximation  for 
the  case  of  a  variable  cross  flow  w(x) . 

2  k 


Demonstration  that 


f  v 

0 


dO  =  0  for  slender  body  section,  - 


The  lift  equation  (lO)  is  particularly  simple  when  the  second  term 

reduces  to  zero.  This  reduction  occurs  when  either  R  -  constant  (the 

2k 

section  under  consideration  is  cylindrical)  or  w  dO  «  0.  It  will 

2k 

f  v 


2it  0 

now  be  shown  that  J  w  d0  =  0  is  obtained  when  the  afterbody  is  suf- 

0 

ficiently  slender  that  slender-body  theory  may  be  applied. 


The  upwash  may  be  obtained  by  inserting  equation  (13)  into  equa¬ 
tion  (8)  and  performing  the  indicated  differentiation.  The  result  is 


w  =  £ua  +  w(xj]  R2/r2(cos20  -  sin20) 
or  on  the  body  (r  =  R), 

w  =  £ua  +  w(x)J  (cos20  -  sin20) 


(15) 


(16) 


which  immediately  leads  to  the  result 

2rt 


f 


w  d0  =  0 


wherever  slender-body  theory  is  applicable. 
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Applications 


Slender  pointed  tody  of  revolution.  -  If,  in  addition  to  being 
slender  everywhere,  the  body  of  revolution  tapers  to  a  point  at  the 

2it 

forward  end,  then  J"*  w  dfl  =  0  and  (p-^  =  0.  Thus  the  lift  equa¬ 


tion  (lO)  reduces  to 


L  = 


2jTqR2<P2 


(17) 


U  sin  9 

With  the  slender-body  value  for  <p9  (equation  (14)),  the  final  result  is 


L  =  2qA2 


(18) 


where  A2  is  the  area  of  the  base  of  the  body.  Thus  the  lift  of  a 
pointed  slender  body  in  an  upwash  field  depends  on  only  the  upwash  and 
the  cross-sectional  area  at  the  base. 


Open- nose  bodies  and  cylindrical  afterbodies  of  bodies.  -  For  a 
circular  cylinder  (dR/dx  «  0),  the  lift  from  equation  (lO)  becomes 


L  _  2qjtR  ^2  ~  *1 
U  sin  9 

If  the  potential  at  the  base  of  the  cylinder  is  given  by  the  slender- 
body  value  (equation  (l4)  with  w  «  0),  the  lift  is  the  same  as  that  of 
a  slender- pointed  body  with  the  same  base  area  because  is  zero. 

Indication  of  the  extent  to  which  use  of  the  slender-body  value  is 
justified  may  be  obtained  from  reference  5  where  an  exact  solution  of 
the  linearized  differential  equation  (equation  (l) )  has  been  obtained 
using  operational  methods  for  the  flow  past  finite  bodies.  Ward’s 
results  for  the  lift  of  a  cylinder  of  length  indicate  that  the 

slender-body  lift  is  realized  within  2  percent  for  Z0/pR>  8. 

The  lift  of  a  cylindrical  afterbody  of  a  body  is  given  by  equa¬ 
tion  (19)  where  is  the  potential  on  the  body  at  the  start  of  the 

aftersection.  The  lift  of  the  afterbody  would  be  zero  if  <p-^  and 
<P2  were  both  taken  to  be  the  slender-body  values  because  would 

then  equal  To  obtain  a  more  accurate  value  for  the  lift  of  the 

afterbody,  <p^_  should  be  obtained  by  more  exact  methods,  such  as  the 
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stepwise  procedure  of  reference  3.  In  the  EXAMPLES  section,  the  lift 
of  a  cone-cylinder  model  is  computed  as  a  function  of  the  cone  angle  by 
this  procedure. 


If  in  the  vicinity  of  the  lip  the  body  is  assumed  to  be  quasi- 
cylindrical  (cLR/dx  «  0),  the  second  term  in  equation  (10)  will  be 
approximately  zero;  if  the  body  is  slender,  the  term  will  still  be  zero 


farther  downstream  because 


f 


w  d0  =  0  as  shown.  Lighthill  points 


out  in  reference  6  that  the  discontinuity  at  the  lip  of  an  open-nose 
body  affects  only  the  region  directly  behind  it  for  a  distance  of  the 
order  of  the  diameter  of  the  body.  The  flow  then  behaves  as  if  the 
open-nose  body  were  a  pointed  body  with  a  continuous  slope.  By  use 
of  operational  methods,  reference  6  shows  that  the  lift  of  an  open-nose 
body  of  sufficiently  large  fineness  ratio  is  given  by  the  slender-body 
value . 


Afterbodies  of  wing-body  combinations.  -  The  analysis  used  for  a 
■body  of  revolution  may  be  considered  to  apply  equally  well  to  a  wing- 
body  combination  provided  the  potential  <p  is  no  longer  restricted 
to  the  form  sin  9  G(x,r).  Consequently  all  the  earlier  formulas  not 
depending  on  this  restriction  still  apply.  In  particular,  equation  (9) 
must  now  be  used  in  place  of  equation  (lO) •  The  lift  of  the  portion 
of  the  body  behind  the  wing  trailing  edge  of  a  wing-body  combination 
may  be  determined  from  equation  (9)  with  <p-^  equal  to  the  potential 

on  the  body  at  the  wing  trailing  edge.  In  order  to  satisfy  the  basic 
assumption  of  constant  interference  downwash  and  sidewash  in  each 
transverse  plane  around  the  body  (due  in  this  case  to  the  wing  and  its 
wake),  the  diameter  of  the  afterbody  should  be  small  compared  with  the 
wing  span.  Reference  7  gives  the  complex  stream  potential  for  a  slender 
body  of  revolution  with  small-aspect- ratio  triangular  wings  mounted  on 
the  cylindrical  part  of  the  body.  The  velocity  potential  is  obtained 
by  evaluating  the  complex  stream  potential  on  the  body  since  the  stream 
function  is  zero  on  the  body.  The  result  is 


^  =  Ua 


(20) 


where  b  is  the  maximum  semispan  of  the  wing  and  R  is  the  radius  of 
the  cylindrical  part  of  the  body.  The  slender-body  approximation  to 
the  surface  potential  at  the  base  of  the  afterbody  is  (equation  (14) ) 


* 


«P2  =  (Ua  +  v2)  sin  ® 
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.where  w„  is  the  induced  upwash  velocity  generated  by  the  wing.  These 

dx 

values  of  <p^  and  <p g  may  now  he  substituted  into  the  first  term  of 
equation  (9)  for  the  lift.  The  second  term  in  equation  (9)  is  zero 

2rt 

•because  r  w  d0  vanishes  as  a  consequence  of  the  slenderness.  The 
0 

result  for  the  lift  of  an  afterbody  of  a  slender  wing-body  combination 
is  found  to  be 

L  -  2qa  [a,  +  A2  (l  +  g).  -  ^!)  -  *(*  ♦  (£)]  (*U 

The  lift  of  a  cylindrical  afterbody  for  the  limiting  case  where 
the  vortex  sheet  is  displaced  sufficiently  far  from  the  afterbody 
that  the  downwash  induced  by  the  wing  at  the  base  of  the  afterbody  is 
negligible  (w2/ba,  »  0)  is  compared  in  figure  2  with  the  lift  of  a 
slender  wing-body  combination  with  no  afterbody.  The  assumption  of 
small  R/b  is  violated  in  the  region  where  the  radius  of  the  body 
approaches  the  semispan.  However,  the  lift  of  the  afterbody  approaches 
zero  in  this  region  and  since  the  lift  given  by  equation  (21)  also  goes 
to  zero,  equation  (21)  may  be  considered  a  good  approximation  to  the 
lift  over  the  complete  range.  The  lift  of  the  wing -body  combination 
is  obtained  from  reference  4  and  is  presented  as  the  ratio  of  the  lift 
of  the  combination  to  the  lift  of  a  pointed  low-aspect -ratio  wing.  The 
lift  of  the  wing  is 

=  2a,jt  (22) 

Combinations  of  bodies.  -  Because  the  linearized  differential 
equation  of  motion  is  assumed  to  apply,  the  solution  for  a  combination 
of  bodies  is  a  linear  superposition  of  the  solutions  for  the  bodies 
alone  and  for  the  interference  effects  between  the  bodies.  The  boundary 
condition  to  be  satisfied  is  that  the  normal  velocity  at  the  surface 
of  the  bodies  is  zero. 

In  order  to  determine  the  interference  lift  for  the  bodies  in 
combination  it  is  necessary  to  know  the  disturbance  one  body  produces 
on  the  other.  In  particular,  the  upwash  and  pressure  field  of  one 
body  contribute  to  the  lift  of  the  other  body.  The  interference  pres¬ 
sure  field  can  be  computed  from  the  interference  upwash  field  by 
integration  of  the  irrotationality  condition 

5u  bw 
bzbx 

to  determine  uj  the  pressure  is  obtained  from  equation  (3). 
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Computation  of  interference  upwash 


Several  methods  may  he  used  to  calculate  the  upwash  field  generated 
by  a  body  at  zero  angle  of  attack  and  at  angles  of  attack.  For  a 
body  at  zero  angle  of  attack  the  method  of  reference  8  may  be  applied 
to  determine  the  upwash.  The  velocity  potential  for  the  supersonic 
flow  about  a  body  at  zero  angle  of  attack  is  given  by 


X 


y (x  -  g)2-32r2 


(23) 


where  f(x)  is  the  local  strength  of  the  source  distribution.  The 
integral  expression  for  the  upwash  is 


w 


sin  6 
r 


f'(Q(x  -  pag 
^(x  - 


(24) 


A  first  approximation  for  the  source  strength  of  a  slender  body  is 


f(x)  =  UR  —  (25) 

dx 

Equation  (24),  with  the  source  strength  given  by  equation  (25),  can  be 
integrated  directly  for  a  given  body. 

An  analogous  procedure  may  also  be  applied  for  nonzero  angle  of 
attack;  that  is,  the  slender-body  value  of  g(x)  (equation  (12)  with 
w  =  0)  for  a  given  body  is  substituted  into  the  integral  expression  for 
the  upwash 

w  _  sin29  H  ^  g’(€)(x  -  g)2dg  +  cos20  P  ^  g(Q(x  -  €)dg 

r2  Jo  Vt*  -  «>z  -  Vtx  -  «2  -  s2--2 

(26) 

and  the  integration  is  performed  without  letting  the  radius  approach 
zero  as  before.  The  result  is  thus  dependent  on  Mach  number,  whereas 
the  upwash  computed  by  the  slender-body  approximation  (equation  (15) 
with  w  =  0)  is  not. 

A  more  accurate  evaluation  of  the  distribution  functions  f(x) 
and  g(x)  is  required  for  bodies  not  necessarily  slender.  For  such 
bodies  a  method  is  presented  in  references  3  and  9  for  evaluating  the 
upwash  by  a  stepwise  procedure,  and  the  stepwise  procedure  is  extended 
to  open-nose  bodies  in  reference  10.  For  a  parabolic  body,  a  comparison 
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is  presented  in  the  EXAMPLES  section  of  the  incremental  upwash  due  to 
angle  of  attack  as  computed  by  slender-body  theory  (equation  (15)  with 
w  =  o),  by  the  more  accurate  stepwise  procedure  of  reference  3,  and  by 
the  intermediate  procedure  of  using  the  slender-body  approximation  for 
g(x)  in  the  integral  expression  for  the  upwash. 

Lift  due  to  interference  upwash 

The  calculation  of  the  portion  of  the  lift  on  a  slender 
body  due  to  the  known  upwash  field  w  from  another  body  is  similar  to 
the  afterbody  problem  in  that  the  upwash  varies  along  the  body.  The 
lift  due  to  w  is  obtained  in  the  same  manner  except  that  <j^  »  0. 

The  upwash  in  each  transverse  plane  is  assumed  constant  and  equal  to 
its  value  at  the  center  line.  The  lift  is  then  given  by  equation  (10) 
with 

<Po  =  Ua  ^  b  sin  G 
Ua  2 


and 


according  to  slender-body  theory.  The  lift  is  thus 


L  =  2qaAo  —  (27) 

Ua 

The  same  result  is  obtained  for  an  open-nose  body  in  a  varying  upwash 
field  if  the  body  is  quasi- cylindrical  (dR/dx  »  0)  near  the  lip. 

The  internal  lift  for  an  open-nose  body  of  length  substantially 
greater  than  p  times  the  diameter  with  supersonic  flow  at  the  inlet 
is  found  from  momentum  consideration  to  be 


L  *  2oqA-j_ 

provided  there  are  no  internal  losses  and  the  flow  is  discharged  at 
the  free-stream  static  pressure.  Generally  this  will  not  be  the  case 
and  terms  must  be  included  to  correct  for  the  internal  axial  forces. 

The  lift  acts  essentially  at  the  nose  and  the  internal  lift  for  an 
open-nose  body  in  a  varying  upwash  field  is  obtained  simply  by  replacing 
the  angle  of  attack  by  the  ratio  of  the  upwash  to  free-stream  velocity 
at  the  nose.  Thus  the  internal  lift  is 
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w-. 

L*  2a-  qA-t  (28) 

Ua  ± 

In  this  analysis  it  is  assumed  that  the  portion  of  the  lift  of  a 
body  that  results  from  the  upwash  generated  by  a  second  body  due  to 
the  presence  of  the  first  body  is  negligible;  in  other  words,  only  the 
primary  interference  effects  are  considered.  In  the  section  EXAMPLES 
the  secondary  interference  effects  are  examined  for  a  cylinder- cylinder 
combination. 

Buoyant  lift  due  to  interference 

The  buoyant  portion  of  the  lift  of  a  body  due  to  the  variable 
pressure  field  from  another  body  is  obtained  from  the  second  integral  B 
of  equation  (5).  Expanding  the  x- perturbation  velocity  u  of  the 
pressure  field  in  a  Taylor fs  series  about  the  horizontal  plane  through 
the  center  line  of  the  body  in  question  yields 


(29) 


where  Zq  is  the  z- coordinate  of  the  center  line  of  the  body.  The 
terms  in  equation  (29)  involving  even  derivatives  do  not  contribute  to 
the  lift,  because  they  are  symmetric  about  the  z  =  zQ  plane.  If  it 

is  assumed  that 


then  that  part  of  the  u-velocity  which  contributes  to  the  lift  is 
approximately 


rVn 

With  the  aid  of  the  irrotationality  condition 
expression  becomes  °z 


u  as  B  sin  6 


rW 

the  previous 
ox 
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Inserting  the  value  of  u  into  the  expression  for  the  buoyant  lift 
(term  B  in  equation  (5))  and  integrating  with  respect  to  6  yield, 
with  w  =  w(x)  , 

X2 

B  =  f  R2 

U  v  \dx 


If  the  body  is  a  cylinder  or  a  quasi  cylinder, 
to  the  pressure  field  is 

B  -  **(?  .  1) 
\Ua  Ua/ 

or 

Wb-K1-1) 


where  A  is  the  mean  cross-sectional  area. 

Equations  (30)  and  (3l)  were  obtained  by  expanding  u  in  a 
Taylor's  series^  hence  they  are  not  valid  in  regions  where  a  discon¬ 
tinuity  in  the  u- velocity  is  present  such  as  occurs,  for  example,  at 
the  intersection  of  the  lip  shock  from  an  open-nose  body  with  another 
body.  The  buoyant  force,  however,  may  be  determined  by  applying  equa¬ 
tion  (30 )  or  (31),  as  the  case  may  be,  a  suitable  distance  behind  the 
intersection*  then  the  buoyant  force  on  the  excluded  area,  associated 
with  the  shock,  is  evaluated  separately  and  added  thereto.  The  buoyant 
force  associated  with  the  shock  is  approximated  herein  by  integrating 
the  pressure  as  if  it  acted  on  the  projected  area  (see  following  sketch) 
of  the  intersection  of  the  shock  with  the  body  on  the  plane  normal  to 
the  plane  passing  through  the  center  lines  of  the  bodies,  and  then 
resolving  the  force  in  the  lift  direction. 


the  buoyant  force  due 

(31a) 

(31b) 


\=z0 


dx 


(30) 
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It  is  assumed  that  u  and  w  are  constant  over  the  shaded  surface 
area  in  the  preceding  sketch  and  have  the  values  appropriate  to  the 
center  line  point  P  immediately  behind  the  shock.  This  assumption 
implies  that  the  body  radius  is  relatively  small  compared  with  the 
separation  of  the  bodies.  For  a  slender  body,  the  projected  area  is 
an  ellipse  with  axes  R  and  PR.  By  virtue  of  the  assumption  on  u, 
the  effective  pressure  acting  on  the  area  is  constant,  having  the  value 
given  at  the  center  line  on  the  downstream  side  of  the  shock  and  zero 
on  the  upstream  side.  The  value  of  the  pressure  on  the  downstream  side 
of  the  shock  is  determined  in  appendix  C  using  the  same  procedure  as 
that  used  in  reference  11  for  the  case  of  the  zero  angle  of  attack. 

The  pressure  coefficient  is 

f  -  -  T M  5in  9°  (32) 

where  the  subscript  0  refers  to  the  open-nose  body.  Since  the  projec¬ 
ted  area  on  which  this  pressure  acts  is  p*®2,  the  buoyant  force  due  to 
the  discontinuity  in  u  at  the  shock  intersection  with  the  body  is 


(33) 
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Equation  (33)  together  with  equation  (30)  or  (3l)  (which  apply  beyond 
the  shaded  region  in  the  sketch)  yields  the  total  buoyant  lift  when  a 
discontinuity  in  the  u- velocity  is  present. 

If  the  body  is  a  cylinder  or  a  quasi  cylinder  in  the  disturbed 
region  of  the  pressure  field  so  that  equation  (31b)  may  be  used,  the 
expression  for  the  buoyant  lift  with  the  lift  due  to  the  wave  dis¬ 
continuity  is  given  by 

,c\  =  As  _  f!i  +  *e  _i_ 

\ L)  B  u  U  q  HR2 

where  A  is  the  projected  area  on  which  the  shock  pressure  acts  and 
Wj  is  the  upwash  immediately  downstream  of  the  wave.  But  for  weak 

waves  the  flow  deflection  and  pressure  coefficient  are  related  by 

Ap  _  2X  _  2 

q  P  P  sin  8q  U 

where  X  is  the  flow  deflection  angle  through  the  wave.  Hence  the 
contribution  of  the  wave  to  the  buoyant  lift  will  be  zero  if 


-1-  =  p  sin  0O  (34) 

itfT 

If  the  effective  projected  area  is  as  assumed  in  the  sketch,  the  rela¬ 
tion  (34)  is  obtained.  Thus,  for  this  case,  the  correct  buoyant  lift 
is  obtained  by  applying  equation  (3l)  with  w^  «  0. 

EXAMPLES 

Lift  of  Cone-Cylinder  Body 

The  lift  of  the  cone  part  of  the  body  as  given  in  reference  3  is 


where  x  is  the  length  of  the  cone.  The  lift  of  the  cylindrical  , 
aftersection  may  be  obtained  from  equation  (19)  with  <p^  equal  to  the 

potential  at  the  base_of  the  cone  and  q>£,  the  slender -body  value 
(equation  (l4)  with  w  =  0} .  The  potential  on  the  body  is 


16 


NACA  TN  2669 


(36) 


Substituting  equation  (14)  with  w  =  0  and  equation  (36)  into  equa¬ 
tion  (19)  and  combining  the  result  with  equation  (35)  yield  the  lift 
of  the  cone-cylinder  which,  expressed  in  coefficient  form  based  on  the 
area  of  the  cylinder,  is 


In  figure  3  the  lift  of  the  cone-cylinder  and  the  lift  of  a  cone 
are  plotted  as  a  function  of  the  parameter  x/pR.  As  would  be  expected, 
the  lift  of  the  cone  and  the  cone- cylinder  approach  the  slender-body 
value  as  the  cone  angle  decreases. 


Interference  Lift  Between  Two  Cylinders  in  Combination 

The  lift  of  two  semi-infinite  cylinders  of  the  same  radius  mounted 
in  the  vertical  plane  and  starting  at  the  same  streamwise  coordinate 
is  computed  to  obtain  some  information  concerning  the  magnitude  and 
the  method  of  computation  of  the  interference  lift  between  bodies. 

There  is  no  loss  in  generality  if  the  cylinders  are  assumed  to  be  in 
the  vertical  plane,  and  the  computation  is  somewhat  shorter.  Also, 
if  a  symmetrically  mounted  vertical  strut  joins  the  bodies,  the  strut 
carries  no  lift  and  produces  no  significant  interference. 

The  lift  of  one  cylinder  in  the  presence  of  the  upwash  field  from 
the  other  cylinder  is  computed  from  equation  (27).  In  figure  4  the 
upwash  at  various  distances  from  a  cylinder  in  the  vertical  plane  as 
computed  by  the  method  of  reference  10  is  presented.  Near  the  lip, 
values  are  improved  by  the  method  given  in  appendix  C.  For  a  short 
cylinder- cylinder  combination,  the  appropriate  value  of  the  upwash  to 
be  used  at  the  base  would  be  the  value  given  by  figure  4.  However,  the 
upwash  approaches  the  slender-body  value  rather  rapidly  downstream  6f 
the  lip  shock.  In  this  example,  therefore,  the  slender-body  upwash 

w2  =  -  UaR2/z02  is  used.  The  lift  due  to  the  upwash  field  in  coef¬ 
ficient  form  and  based  on  the  cross-sectional  area  of  the  cylinder  is 
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“  ~  2a  2 
z0 

where  Zq  is  the  distance  between  the  center  lines  of  the  cylinders. 

The  interference  lift  given  by  equation  (38)  (which  may  be  regarded 
as  the  first  iteration  of  the  total  interference  lift)  can  be  considered 
to  result  from  a  certain  distribution  of  doublets  along  the  axis  of  the 
body.  This  doublet  distribution  produces  an  upwash  field  as  though 
(in  slender-body  approximation)  the  upper  body  were  moving  in  the  nega¬ 
tive  z-direction  with  velocity  aXJ  zq2/h2.  The  upwash  induced  at  the 
lower  body  due  to  this  motion  at  a  distance  Zq  below  the  center  line 
of  the  upper  body  is  a  fraction  z q2/r2  of  this.  Thus  the  lift  of  the 

lower  cylinder  in  the  presence  of  the  upwash  from  the  upper  cylinder 
due  to  the  presence  of  the  lower  cylinder  is  given  by 


(39) 


The  complete  solution  for  the  lift  due  to  the  upwash  for  both  cylinders 
is  then 


The  validity  of  this  equation  deteriorates  as  zQ  diminishes  because 

of  the  neglect  of  the  gradient  in  upwash  across  each  body.  In  the 
limiting  case  when  the  cylinders  are  touching  (r/z^  ~  1 /2),  equa¬ 
tion  (40 )  is  certainly  invalid  quantitatively,  but  the  behavior  of  the 
iterations  is  probably  qualitatively  correct.  For  this  limiting  case 
equation  (40)  becomes 


Cl  =» 


4 

—  a 

5 


(41) 
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The  magnitude  of  the  terms  of  higher  order  than  R2/zQ2  for  the  contig¬ 
uous  case  is  5  percent  of  the  external  lift  of  the  cylinder  alone.  An 
exact  value  for  the  upwash  interference  lift  at  small  separation  Zq,  within 

the  framework  of  slender-body  theory,  may  he  found  by  considering  the 
two-dimensional  flow  around  two  cylinders.  This  flow  is  equivalent  to 
a  doublet  external  to  a  cylinder.  By  use  of  the  method  of  images 
(reference  12),  the  potential  in  series  form  (because  each  reflected 
doublet  violates  the  boundary  condition)  may  be  obtained  and  hence 
the. lift.  The  result  for  the  interference  lift  is 


(42) 


For  the  contiguous  case  the  interference  lift  due  to  the  upwash  field  is 

+  .  .  +  (-l)n+1  - - - -  -  .  .  .1  (43) 

[4  9  16  25  (n  +  l)2  J 

4a(l  -  Jt2/l2) 


=  -  4a 


The  magnitude  of  higher-order  terms  for  this  method  is  about  7  percent 
of  the  external  lift  of  the  cylinder  alone.  The  difference  between 
equations  (43)  and  (4l)  is  due  to  the  fact  that  in  obtaining  equa¬ 
tion  (40)  the  boundary  condition  is  only  approximately  satisfied. 

A  procedure  analogous  to  that  used  to  obtain  equations  (40)  and 
(42)  may  be  applied  to  obtain  expressions  for  the  buoyant  lift.  The 
buoyant  lift  is  equal  to  the  interference  lift  due  to  upwash.  For  this 
example,  the  internal  interference  lift  is  zero  because  the  noses  of 
the  cylinders  are  ahead  of  the  disturbance  fields . 


The  total  lift  of  the  combination,  including  the  internal  lift 
for  the  isolated  cylinders,  with  terms  of  higher  order  than  R2/z^2 
neglected,  is 


=  8a  I  1 


(-a 


% 


(44) 
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The  interference  lift  may  be  expressed  in  terms  of  an  interference  ratio 


i  = 


dCL( 

da 


dC 


da 


da 


da 


dC- 


_^b 

da 


v 


(45) 


where  CT  is  the  lift  of  the  isolated  bodies. 

•4) 

It  is  interesting  to  note  that  the  interference  lift  ratio  for  the 
case  of  the  cylinders  in  the  horizontal  plane  would  be 


i  = 


(46) 


*0 


Interference  Lift  Between  Two  Parabolic  Bodies  in  Combination 

In  this  example ,  the  interference  lift  between  two  identical 
parallel  parabolic  bodies  mounted  in  the  vertical  plane  and  at  incidence 
angle  a  with  respect  to  the  free-stream  is  computed. 

The  upwash  distribution  about  a  parabolic  body  due  to  angle  of 
attack  may  be  computed  either  by  slender-body  theory,  by  the  stepwise 
procedure  of  reference  3,  or  by  the  intermediate  procedure  of  using  the 
slender-body  approximation  for  the  g(x)  in  the  integral  expression 
for  the  upwash  (equation  (26)).  The  equation  of  a  parabolic  body  in 
dimensionless  form  is 

R*  =  \  (l  -  x*)x*  (47) 


where 
R*  PR/Z 

F*  fineness  ratio/p  =  F/p  =  t^^max 
x*  x/Z 

length  of  body  pointed  at  both  ends 


l 


20 


NACA  TN  2669 


The  up-wash  in  the  vertical  plane  through  the  body  center  line  obtained 
by  the  intermediate  procedure  for  the  parabolic  body  defined  by  equa¬ 
tion  (47)  is 


w 

Ua 


where  Zq*  =  PZq/7.  This  equation  is  applicable  between  the  Mach  cones 
from  the  pointed  nose  and  end  of  the  body;  that  is, 


z0*$  X  $  1  +  ZQ* 


The  upwash  as  computed  by  slender-body  theory  is 


w 

Ua 


(49) 


The  comparison  of  equations  (48)  and  (49)  with  the  upwash  computed 
using  reference  3  for  F*  =  7.5  and  12.5  and  z^*  =  0.075  and  0.125, 
respectively,  is  presented  in  figure  5.  It  is  seen  that  the  upwash 
given  by  equation  (48)  varies  in  the  same  manner  with  respect  to  magni¬ 
tude  and  position  as  that  obtained  using  the  method  of  reference  3. 

Since  the  computation  with  equation  (48)  is  faster,  equation  (48)  is 
used  herein  for  the  interference  lift  calculations  due  to  the  upwash 
from  the  interfering  body.  The  sign  of  the  interfering  upwash  acting 
at  the  base  of  the  primary  body  is  important  because  it  determines 
whether  or  not  the  incremental  lift  due  to  interference  upwash  is 
favorable.  The  upwash  given  by  slender-body  theory  will  always  have 
the  same  sign  for  a  given  plane  of  symmetry,  whereas  the  upwash  given  by 
both  of  the  more  exact  methods  changes  sign  (see  fig.  5).  For  the 
calculation  of  the  buoyant  lift,  the  simpler  equation  (49)  (rather 
than  equation  (48))  has  been  used  in  equation  (30)  because,  for  inte¬ 
grated  values,  the  difference  between  the  various  curves  is  unimportant. 


In  figure  6  a  contour  plot  of  the  interference  ratio  is  presented 
for  the  case  where  the  distance  between  the  center  lines  of  the  bodies 
varies  from  z q*  =  0.075  to  z*  =  0.150  and  one  body  moves  downstream 

from  x*  =s  0  to  x*  «  0.3  with  respect  to  the  other  parabolic  body. 
The  length  of  each  body  is  0.87.  The  maximum  lift  curve  slope  occurs 
when  the  body  axes  are  farthest  apart  and  one  body  is  at  the  most 
rearward  station  of  the  positions  considered. 
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Interference  Lift  Between  Open-Nose  Body  and  a 
Parabolic  Body  in  Combination 

A  contour  plot  of  the  interference  ratio  at  a  Mach  number  of  2  is 
presented  in  figure  7(a)  for  an  open-nose  body  mounted  above  or  below 
and  downstream  of  a  parabolic  body.  The  bodies  are  parallel  and  at 
incidence  angle  a  with  respect  to  the  free  stream.  All  pertinent 
dimensions  are  shown  in  the  figure.  In  figure  7(b)  another  open-nose 
body  identical  to  the  first  is  symmetrically  added  to  the  configuration 
and  the  interference  ratio  is  presented  for  this  arrangement. 

For  these  plots,  the  upwash  distribution  about  the  parabolic  body 
was  computed  from  equation  (48)  with  F  equal  to  15.  For  the  open-nose 
body  the  upwash  was  determined  using  the  method  of  reference  10.  For 
the  computation  of  the  buoyant  force,  equation  (3l)  was  used  with  the 
mean  radius  determined  by  the  amount  of  volume  of  the  body  in  the  dis¬ 
turbed  field  at  a  particular  position.  The  lift  of  the  isolated  com¬ 
ponents  was  obtained  using  slender-body  theory  (equation  (l8)  with 
w  *  0). 


Figure  7(a)  indicates  that  for  the  positions  considered,  the 
maximum  lift-curve  slope  occurs  when  the  open-nose  body  is  at  the 
farmost  aft  and  outboard  station  from  the  parabolic  body.  The  same 
result  is  noted  in  figure  7(b)  for  the  parabolic  body  -  two  open-nose 
body  configuration. 

As  a  matter  of  general  interest,  the  interference  lift  at  zero 
angle  of  attack  is  presented  in  figure  8  for  the  one  open-nose  con¬ 
figuration.  When  the  open-nose  body  is  mounted  above  the  parabolic 
body,  the  lift  is  negative;  and  when  it  is  mounted  below,  the  lift  has 
the  same  magnitude  but  the  opposite  sign.  In  the  calculation  of  the 
lift  the  upwash  about  the  parabolic  body  was  determined  from  reference  8, 
and  the  upwash  about  the  open-nose  body  was  determined  from  reference  10. 


Lewis  Flight  Propulsion  Laboratory 
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Cleveland,  Ohio,  December  17,  1951 
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APPENDIX  A 
SYMBOLS 

The  following  symbols  are  used  in  this  report: 

2 

cross-sectional  area  of  Body  of  revolution,  itR 
maximum  semispan 

buoyant  lift  (second  integral  in  equation  (5)) 
lift  coefficient,  lift/qA2 

fineness  ratio  of  parabolic  body  of  revolution, 
length/maximum  thickness 

f unctions  of  x  and  r  (defined  in  appendix  B) 
local  strength  of  source  distribution 
derivative  of  local  strength  of  doublet  distribution 
interference  lift  ratio 

lift  due  to  angle  of  attack  and  upwash  (first  integral  in 
equation  (5)) 

length  of  body  pointed  at  both  ends 

Mach  number 

integer 

deviation  of  local  pressure  from  free-stream  pressure 

principal  part  of  Ap  contributing  to  lift 

1  2 

free-stream  dynamic  pressure,  —  pU 

2 

radius  of  body 
maximum  radius  of  body 
radial  coordiante, 
wing  area 
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t 

U 


u,w 


u,v,w 


integration  variable 
free- stream  velocity 

perturbation  velocities  in  x-  and  z-directions ,  respectively 
perturbation  velocity  components  due  to  interfering  body 


V 

a 

P 

r 

e 

l 

p 

q> 

x 

* 


resultant  velocity 

components  of  resultant  velocity  in  r~,  0-,  and  X-directions, 
respectively 


Cartesian  coordiantes 
angle  of  attack 

cotangent  of  Mach  angle  of  free  stream, 

adiabatic  exponent 

tangential  coordiante,  sin"-^(z/r) 


M2  -  1 


integration  variable  (x-coordinate  of  sources  or  doublets) 
density 

total  perturbation  potential  for  primary  body 
portion  of  <$  contributing  a  nonzero  value  of  lift 
portion  of  #  corresponding  to  zero  angle  of  attack 
portion  of  #  due  to  sidewash 


Subscripts: 
0 


coordinate  measured  from  axis  of  body  generating  disturbance 
field 

start  of  body 
base  of  body 
bodies  alone 
bodies  in  combination 
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i  due  to  interference 

w  wing 

Subscript  coordinates  indicate  partial  differentiation  with  respect 
to  subscript  variable. 

Primes  denote  ordinary  differentiation. 
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APPENDIX  B 

FORM  OF  VELOCITY  POTENTIAL  FOR  BODY  OF  REVOLUTION  AND  PRINCIPAL 

PART  OF  PRESSURE  CONTRIBUTING  TO  LIFT 

The  total  perturbation  potential  for  the  primary  body  is  assumed 
to  satisfy  the  governing  equation  for  small-disturbance,  isentropic 
irrotational  flow 


3^a  -a  -  a  as  o 

H  xx  wyy  ^zz  u 

The  general  solution  of  equation  (biJ  as  given  in  reference  13  is 

ao  OP 

*C*u,e)  =  E  cos  m  9  Fm(x,r)  +  sin  m  6  G^fx^r) 


(Bl) 


m=sO 


UfrsO 


where 


Fm= 


ftir/ 


fffi(x  -  (3r  cosh  t)dt 


cosh-1 

pr 


*  -  r111 


x-pr 

e  in 


fmCOd£ 


0  /^(x  -  £)2  -  p2r2 


x  m  0 

(-  T-j  P  6m(x  "  Pr  costl  t)dt 

',r  TJ  ”  v-1  x 

cosh  — 

Pr 


=  -  i* 


m  X-Pr 

H)  f 


0  ^(x  -  £)2  -  p2r2  J 


}  (B2) 


The  linearized  boundary  condition  on  the  body  is 
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This  boundary  condition  is  correct  to  the  order  a  dR/dx  or  (dR/dx)2. 
Because  of  the  form  of  the  boundary  condition,  the  total  potential 
must  be  given  by 


$  m  sin  0  G^+  cos  0  +  Fq 


“  <P  +  +  +  X 


(B4) 


(with  v  «  v(x),  w  a  w(x)).  The  part  of  the  pressure  contributing  to 
the  lift  will  be  obtained  to  the  order  of  accuracy  of  the  linear  partial 
differential  equation  and  the  linearized  boundary  condition.  The 
exact  form  of  the  pressure  coefficient  is  given  by  the  compressible 
Bernoulli  equation 


AP 

q 


M2 


(■  ■  S)] 


r 

r- 1 


(B5) 


Expanding  equation  (B5)  yields  the  following  approximate  expression 
for  the  pressure  coefficient: 


(B6) 


This  equation  is  correct  to  the  second  order  in  1  -  V2/U2.  The  net 
velocity  at  any  point  is 


where 


V2  *  vx2  +  vR2  +  ve2 

vx  -  u  +  u  +  $x  (B7a) 

Vr  *  (Ua  +  "wjsin  0  +  v  cos  0  +  $r  (B7b) 

—  _  $  Q 

V0=*  (Ua  +  w)cos  0  +  v  sin  6  +  —  (B7c) 


Substituting  equation  (B3)  into  (B7b)  yields  the  net  velocity  in  the 
radial  direction  of  the  body 


vR  =  u  — 


dx 


(B8) 
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Substituting  equations  (B8),  (B7a),  and  (B7c)  into  the  first  term  in 
equation  (B6)  yields,  on  the  body. 


<1 


.  -  f  +  ^  +  5'(1  +  “'\+!j!  +  /™\2  +  /a  +  » 

_U  u2  v  \  U/  jjZ  \dx/  \  u 


7 


cos20  + 


sin  6 


(B9) 


The  following  terms  in  the  approximate  pressure  coefficient  yield 
nonzero  values  for  lift  when  equation  (B4)  is  substituted  into  (B9), 
the  indicated  differentiation  is  performed,  the  result  is  substituted 
into  the  lift  integral 


Lift  = 


2jc 


f 


—  sin  9  R  d6  di 

q 


(Bio) 


and  the  integration  with  respect  to  9  is  performed: 

Sp  _  2<*x  2^xu  2<^x%x  2u  u2 

q  U  U2  U2  u  U2 


(Bll) 


All  other  terms  integrate  to  zero .  If  the  next  term  in  the  expansion 
of  the  pressure  coefficient  (equation  (B6))  is  included,  by  the  same 
procedure 


2<px  t  2p2$xq  t  Z&xXx 

u  +  u2  +  u2 


2u  P2u2 
U  +  U2 


(B12) 


where  the  terms  of  higher  order  than  those  appearing  in  equation  (Bll) 
have  been  neglected.  It  can  be  shown  from  slender-body  theory  that 
*xM  S/Cf,  and  Xx/u  are  of  order  a  dR/dx  or  (dR/dx)2. 


Neglecting  terms  of  higher  order  than  a  dR/dx  or  (dR/dx)2  in 
equation  (B12),  where  a  and  dR/dx  are  considered  to  be  small  and 
of  the  same  order,  yields 


5p  _  2<*x  2u  _  2u  2u 
q  U  U  U  U 


(B13) 
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The  component  potentials  ^  and  X  thus  contrihute  only  higher- 
order  terms  to  the  lift.  _The  surviving  component_in  equation  (B4)  is  <p, 
■which  is  due  to  a  and  w  alone;  the  side-wash  v  does  not  enter. 

Comparison  of  equations  (B12)  and  (B13)  indicates  that  at  high 
Mach  numbers  some  of  the  neglected  terms  may  become  as  large  as  those 
retained. 
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APPENDIX  C 


PERTURBATION  VELOCITY  COMPONENTS  IN  VICINITY  OF  LIP 


SHOCK  FROM  AN  INCLINED  OPEN-NOSE  BODY 


The  method  used  herein  to  estimate  the  velocity  components  in  the 
vicinity  of  the  lip  shock  from  an  open-nose  body  at  angle  of  attack  is 
'  similar  to  the  technique  used  in  reference  11  for  the  linearized 
treatment  of  the  lip  shock  from  an  open-nose  tody  at  zero  angle  of  attack. 
The  boundary  condition  on  the  open-nose  body  in  integral  form  may  be 


written 


x-pR  _ 

£i2_£  P  g' (*>)(*  ~  gLgL-  *  .  Ua  sin  e 

r2  o  V(x  -  s)2  -  p2r2 


(Cl) 


where  the  lip  is  located  at  x  a  fR, 
x  -  pR  is  very  small  compared  with 
to 


UcxVHr 

o2 


-r 


In  the  vicinity  of  the  lip, 
pR  and  thus  equation  (Cl)  reduces 


x-PR 


g'(€)d.g 


*Jx  -  PR  -  £ 


(C2) 


This  is  in  the  form  of  Abel's  integral  equation  for  the  unknown  function 
g'(x)$  the  solution,  according  to  reference  14,  is 


g'(x) 


(C3) 


when  R  is  assumed  to  be  approximately  constant  in  the  vicinity  of  the 
lip.  Introducing  equation  (C3)  into  the  expression  for  the  upwash 
(equation  (26)),  restricting  consideration  to  the  vicinity  of  lip  shock 
(x  -  pr  very  much  less  than  pr) ,  and  integrating  yield 

JL  =  -  sin20  +  cos2©  0L  -  l)  (C4) 

In  the  integration  of  the  last  term  in  equation  (26),  the  function  g(x) 
is  given  by  the  integral  of  equation  (C3),  which  is 

g«  -  yj  +  C  (05) 

P2rt 

with  the  constant  C  equal  to  zero  because  g(0)  =*  0 . 
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The  u-velocity  component  in  the  vicinity  of  the  shock  is  obtained 
by  inserting  equation  (C3)  into  the  integral  expression  for  the  u- 
velocity  and  neglecting  x  -  |3r  compared  with  pr  as  before.  The  result  is 

.  ^  ^  ^  sin  0  (C6) 

U  q  P  Vr  v 

At  the  lip  shock  the  upwash  is  zero  when  calculated  by  the  method  of 
reference  10.  In  figure  3  the  upwash  given  by  equation  (C4)  is  joined 
by  a  dashed  line  to  the  values  given  by  the  method  of  reference  10 
downstream  of  the  lip  shock.  The  slope  of  the  dashed  line  is  about  the 
same  as  that  given  by  reference  10,  thus  indicating  the  validity  of  the 
method  of  this  appendix. 
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>  Figure  2.  -  Lift-curve  slope  ratio _for  wing-tody  -  aftertody 
combination  with  Wg/Ua  =  0, 
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Figure  3.  -  Comparison  of  lift  of  cone-cylinder  with  lift 

of  cone. 
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isolated 

components 


Figure  8,  -  Contour  plot  of  interference  lift  coefficient  at  zero  angle 
of  attack  for  open-nose  body  mounted  above  parabolic  body. 
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afterbody  of  a  slender  wing-body  combination,  and  afterbody  of  a  slender  wing -body  combination,  and 

three  combinations  of  bodies  to  illustrate  the  method.  three  combinations  of  bodies  to  illustrate  the  method. 


